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Abstract: We study the phenomenology of, and search techniques for, a class of "Hidden 
Valleys." These models are characterized by low mass (well below a TeV) bound states 
resulting from a confining gauge interaction in a hidden sector; the states include a spin-one 
resonance that can decay to lepton pairs. Assuming that the hidden sector communicates 
to the Standard Model (SM) through TeV suppressed operators, taking into account the 
constraint from the Z pole physics at LEP, searches at Tevatron may be difficult in the 
particular class of Hidden Valleys we consider, so that we concentrate on the searches at 
the LHC. Hidden Valley events are characterized by high multiplicities of jets and leptons 
in the final state. Depending on the scale of confinement in the hidden sector, the events 
are typically more spherical, with lower thrust and higher incidences of isolated leptons, 
than those from the SM background processes. Most notably, high cluster invariant mass 
and very narrow, low mass resonances in lepton pairs are the key observables to identify the 
signal. We use these characteristics to develop a set of cuts to separate the Hidden Valley 
from SM, and show that with these cuts LHC has a significant reach in the parameter 
space. Our strategies are quite general and should apply well beyond the particular class 
of models studied here. 
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1. Introduction 

In the anticipation of major discovery at the Large Hadron Collider (LHC) for physics 
beyond the standard model (SM) at the TeV scale, it is prudent to keep in mind that the 
high energy frontier may open up the possibility for new physics with non-conventional 
forms. Many theoretical models include a "hidden sector" that couples to the SM weakly. 
The hidden sector may communicate to the standard model through a heavy mediator 
(typically with a mass scale of TeV or higher), charged under both the SM and hidden 
sector. TeV scale particle physics model building is replete with examples of such hid- 
den sectors. These include gravity mediated SUSY breaking [1-3], gauge mediated SUSY 
breaking [4-6], Twin Higgs models [7], milli-charged hidden sector dark matter [8,9] or 
unparticle physics [10,11]. String models also generate hidden sector matter [12,13]. 

We consider in this paper a class of hidden sector models which are characterized by 
a new confining gauge interaction in the hidden sector. The dynamics in the hidden sector 
is set by the confinement scale, which introduces a mass gap into the theory. The addition 
of new quarks into the hidden sector will give rise to a variety of hidden bound states (as 
in QCD), with presumably low mass much below 1 TeV. This class of models, introduced 
in [14], in an example of a "confining Hidden Valley." The term "Hidden Valley" refers 
to the presence of the low mass "valley" (v-)states (f-quarks and iJ-hadrons in the present 
case), which can only be observed by passing through or over a barrier separating the 
hidden sector from the standard model sector. The novel phenomenology of the Hidden 
Valley arises because of the low mass of the u-particles and the nontrivial dynamics in the 
hidden sector. The masses of the f-hadrons are determined by either the u-quark mass 
or the confinement scale, whichever is larger. The lower bound on the ii-hadron's mass is 
derived by requiring that it decay before BBN; the precise constraint thus depends on the 
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dimension of the operator mediating communication between the hidden sector and the 
standard model, but it is typically sub-GeV. When the mass of the bound state is below 
10 MeV, on the other hand, the dominant constraints may be derived from astrophysics 
and cosmology. 

Similar to the mediation of SUSY breaking, the effects of a Hidden Valley on our SM 
observable sector are crucially determined by the nature of the mediators [14-16]. Given 
the presumably large gap between the mass of t;-scctor particles and the heavy mediators, 
there is effectively a barrier (like a mountain pass) between the SM and Hidden Valley 
which must be surmounted in order to produce the valley particles via the SM processes. 
When the mediator is integrated out, an effective operator is generically of the form 

OyOsM 

9v9SM jf^k ' ^^-^^ 

where gy and ggj^ are the mediator couplings to the two sectors, M the mass of the 
mediator and the power of the suppression, k, depends on the mediator and the nature 
of the hidden sector. Examples of possible mediators are Z"s, Higgs sectors, messenger 
fermions which connect the gluons from the hidden and standard model sectors, or gravitons 
from Randall-Sundrum or large extra dimensions. On the other hand, if the mediator scale 
M is accessible by the experiments, one could thus expect enhanced effects. 

Once the v-sector quarks are produced through the mediator, they confine into v- 
hadrons. If the f-hadrons arc much lighter than the mediator, many of them are produced 
in the v-sector hadronization. Provided the mass gap is not too low and the hidden sector 
resonances not too light, these v-hadrons decay back through the mediator to combinations 
of leptons and quarks which, when summed together, are standard model neutral since the 
i;-hadron itself carries no standard model charges. Thus i;-hadrons, including the ones that 
will be the focus of our study, may decay to neutral pairs of leptons, heavy quarks, or light 
quark jets. In particular, the decay of spin-one resonances tends to be democratic among 
standard model fermions. 

To be concrete and to gain a basic understanding of the features and search strategies 
at a collider, we choose a particular model for the hidden sector and mediator. We take a 
Z' mediator and a hidden sector with one light quark, which has narrow spin-zero and spin- 
one resonances and no absolutely stable ?;-hadrons. Another example of such a model has 
been studied elsewhere [17] though no signal-to-background studies have been performed. 
Many of the features we describe here will be transferrable to a broader class of models 
with a confining gauge group in the hidden sector. We choose to focus on the region of 
parameter space in this model where the v-hadrons decay promptly back to SM particles, 
which was shown in [15] to occur for 

m^^>30 GeV. (1.2) 

Lighter mass i;-hadrons will either decay with a displaced vertex in the detector or will 
result in missing energy, leading to qualitatively different signals which are not the focus 
of this study. The general features we will study for the prompt decays are independent, 
in a qualitative way, of the mass of the v-hadron and the specific nature of the mediators. 
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Table 1: Charges of the SM and hidden sector quarks under the mediator gauge group, C/(l)^, 
and the hidden sector confining group, SU{Nc). The hidden sector quarks are uncharged under 

SU{3) X SU{2) X [/(f) and q+ + ^ -2. 



We will demonstrate that these new confinement dynamics may lead to quite distinctive 
features that are characteristically different from the SM expectations, and demonstrate 
the ways in which the hidden sector strong dynamics may be separated from SM QCD. 

It should be kept in mind that Hidden Valley phenomenology can be diverse far beyond 
the simple setting laid out above. It depends on the size of the mass gap, the nature of 
the mediator, and the matter content of the hidden sector, all of which enter the effective 



operator, Eq. (1.1). For example, in the limit that the mass gap is taken to zero and the 
quarks are massless, we recover a model of scale invariance, similar to Unparticles [10,11]. 
On the other hand, if all the quarks are much heavier than the confinement scale, v- 
quarkonium states can play an important role in the physics [14], and if some of the heavy 
t;-quarks also carry standard model charges, "quirk" phenomenology results [18, 19]. In 
some models, certain of the bound states may be stable leading to large missing energy 
signals. Some of the states may be quasi long-lived, decaying a macroscopic distance away 
from the interaction region giving rise to the displaced vertex discussed in [14,15]. With 
a new Higgs mediator that mixes with the SM Higgs boson, Hidden Valleys may give rise 
to novel search techniques for the Higgs through highly unusual decay patterns, some of 
which may be implemented at the Tevatron and LHC. Clearly there is a rich and broad 
phenomenology which remains to be explored. 

The rest of paper is organized as follows. In Sec. 2, we outline a general scenario of 
a Hidden Valley model, reiterate the current bounds, and calculate the typical production 
rates for the signal at the LHC. Our treatment for the w-quark hadronization is also de- 
scribed. In Sec. 3, we systematically explore a method for separating Hidden Valley signals 
from the QCD background which relies on distinguishing the unique Hidden Valley event 
shapes; we discuss the predominant backgrounds to our signal. We conclude in Sec. 4. 



2. The Model And The Hidden Valley Particle Production 
2.1 Model description 

We utilize the model of [14] in computing production and decay. In this model, SM 
fermions such as q and q annihilate through the Z' mediator to hidden sector quarks v and 
V. The charges are repeated for convenience in Table |l[ We choose the simplest hidden 
sector content possible, a single light u-quark, which we denote U in the table, though an 
additional heavy hidden sector quark (not relevant for the phenomenology discussed here) 
must be added to make the theory anomaly free. 
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Figure 1: (a) Production of w-quarks and (b) decay of ?;-hadrons, pv (the vector bound state of 
w-quarks), rj^ (the pseduoscalar bound state). 



Since the Z and Z' are charged under both the hidden and visible Higgs sectors, the 
Z and Z' mix through mass terms. As a result, t;-quarks may be produced through either 
a. Z OY Z' mediator, as indicated in Fig. |l|a. The most stringent constraints on models of 
this type arise from the LEP measurements on the Z pole. We assume that no more than 
a few events of u-quark production will be consistent with the constraints from LEP, so 
that we require a branching fraction less than 5 x 10^'' from decays of Z to w-quarks at 
LEP through Z — Z' mixing. In particular, the production cross-section is 



4!^V^Vl + 2!^L (2.1) 



{s — m^^iY ^V\,rr?^, \ s 

where g' is the gauge coupling of the Z', is the charge of the u-quark under U{l)y., 
is the mass of the u-quark, and Rf includes the effects oi Z — Z' mass mixing: 

nL,R ^L,R z' r,r,L,R Zn Co 9^ 

•' s — + iL z^z 

This is to be compared against the standard model total cross-section: e^e" ^ Z — > 
anything, asM = dzf^z / z- Requiring a^jasu < 5 x 10^^ imposes a bound niz'/g' ^ 
7 TeV. In the remainder of this paper, we take 

mz' = l TeV, mz'/g' = 7 TeV, iV^ = 3. (2.3) 

Later we generalize the parameters to determine the LHC reach. 

2.2 On-threshold production: The v-onia production 

Historically, new bound-states have been discovered on threshold, where ^/s = 2mq. The 
existence of c and b quarks was gleaned through the observation of the J/ ip and T bound- 
states. It is thus natural to consider the vv bound state formation as the signal production 
mechanism. 
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Figure 2: Total cross section of u-liadron production versus the w-quark mass m„ at (a) tlie 
Tevatron and (b) tlie LHC. Tlie peaked curves are on-tlireshold vv bound-state production, the 
flatter curves open-flavor vv pair production. Here and throughout this paper the parameters 
niz' — 1 TeV, niz' / g' — 7 TeV and iV^ = 3 arc chosen. 



For the case of a heavy quark, the calculation of on-threshold bound-state production 
follows that of quarkonium (see [20] for details of quarkonium production): 



f^J'*7.[\R^\\\RR\'){Ql+Q%)\4>^^,\'m,n 



(rthresh= I ^ dxi I ^ dX2Tf^ (Sx.xo-J .)^+Vl.ml. KSxiX2 - m^^) 



xE,2/,(xi)/5(x2), (2.4) 

where S is the hadronic center of mass energy (5x1X2 = s). For perturbative quarkonium, 
the wavefunction cp^fi can be computed analytically in a manner analogous to the hydrogen 

1 3 

^N~^mvas{ml) /it, where my is the f-quark mass. For non- 



wavefunction: |</>^/i|^ 



perturbative bound states of vv, the wave function cannot be computed analytically, but 
can only be approximated through the relation |i<V_Ap = 8^c'rnvh\4'vh\'^ , with Fy^A the 
vector and pseudoscalar decay constants. The decay constants must be determined either 
experimentally or by non-perturbative lattice calculations. 

The production cross-section is shown in Fig. ^ as a function of the t;-quark mass (a) 
at the Tevatron and (b) at the LHC energies for a mediator mass niz' = 1 TeV, assuming 



- 5 - 




Figure 3: Differential cross-sections as a function of partonic cm. energy y^. The solid curve 
represents the inclusive cross-section without acceptance cuts, the dashed is with the trigger re- 
quirement in Eq. (3.3), and the dotted includes further cuts as given in Eqs. (3.7), (3.10), and 

(EH- 



Fv,A = "i^/i (the result can be easily rescaled for different decay constants). We see the 
large resonant enhancement at 2mv ~ mz or mz' due to the mediators between the SM 
fields and the f-sector. Off the Z or Z' resonances, the production cross-sections are 
typically of the size a few fb at the Tevatron and lO's of fb or smaller at the LHC. The 
smallness of the cross section is largely due to the suppression of the Z' propagator. Given 
the several orders of magnitude of larger rates for the SM DY production as backgrounds, 
the f-onia signal may not be readily observable. We will not pursue this channel further. 

2.3 Open-flavor production and hadronization 

Once produced, the v-quarks will hadronize into a "f-color" singlet state, as in QCD. If 
they are produced near the threshold 2m^, it is likely that the vv will form a single bound 
state, pv or rjv (the vector and pseduoscalar bound states of w-quarks, respectively), as 
we discussed in the previous section with production rates on-threshold shown in Fig. |2|. 
If produced well above threshold, they form multiple u-color singlet bound-states, leading 



to possible f-jets. We make use of the production cross-section Eq. (|2.lD , appropriately 
modified with parton distribution functions at a hadron collider. The total open-flavor 
cross-section is shown in Fig. |2| as a function of (the flatter curves). As expected, 
the production rate is much larger than the bound state production when 2my ^ mz or 
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mz' ■ The u-quarks are now being produced through a resonant decay of the mediator that 
couples the qq in the initial state to vv in the final state. The cross section is nearly constant 
up to the kinematical limit rriy ~ mz/'^, razi/2, and the step-function nature reflects the 
fact crossing over the mediator thresholds. Once again, we see that the contribution from 
a 1 TeV Z' is significantly larger than that from the Z at the LHC, while it is opposite 
at the Tevatron due to the limited parton cm. energy. The cross section scales with the 
model parameters as 

{g'"^ /mz' at the Z' resonance, 

{g' /niz')'^ off the ^'resonance, even at the Z resonance. 

Given the expected integrated luminosity of a few fb^^ at the Tevatron, the signal rate 
would be rather low, and it will be thus a challenge to observe the Hidden Valley effects. 
We will henceforth focus on the search at the LHC. 

At the LHC, the majority of the contribution comes from Z', three times as much 
as that from Z. The characteristic kinematics to look for the vv is at the invariant mass 
peak of the sub-process, as shown in Fig. |3|, where a 1 TeV Z' and mixing with the SM 
Z are evident. However, whatever the mediator turns out to be, it is a character with 
known properties. We thus will make use of the knowledge about it such as the mass, 
width and typical couplings to the SM particle etc. to deduce the existence of the Hidden 
Valley signal. To simplify the analyses, we will devise an initial cut to center on the new 
mediator mass, as most of the rate is produced there. For definiteness and convenience for 
the presentation, we choose 

mz - 50 GeV < < m,^ + 50 GeV for Z mediator (2.6) 
mz' - 100 GeV <y/s<mz' + 100 GeV for Z' mediator. (2.7) 

This pre-selection condition will effectively separate the background considerations as well. 
For the Z', this requirement is essentially a lower cut on the energetics, given the steeply 
falling parton luminosities. 

To incorporate the hadronization process for the pair-produced v quarks, we write 



(^vh = j dxidx2da{s)qq^vv Pvv^h's, (2. 



where da{s)qq^vv is differential cross section for the partonic process qq vv, and Py^^h's 
is a production probability for the initial 7;u-quark pair to transform into w-hadrons, which 
is described by hadronization models, such as Lund Fragmentation Model [21], Webber 
Cluster Model [22], and Quark Combination Model (QCM) [23-27] etc. Here we adopt a 
modified QCM where the baryon production is neglected. Within this model, A^-pairs of 
new vv can be produced from the vacuum just as standard model light quark pairs. N is 
assumed to satisfy a Poisson Distribution: 

P((Ar),iV_i) = ^_^e-W (2.9) 
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Figure 4: Probability distribution of an event containing N w-hadrons (a) via a Z mediator and 
(b) via a Z' mediator. 

where (A^) is the average number of quark pah's. According to Ref. [23,25], (A^) is param- 
eterized by 

(N) = \Ja^ + 13^ -a - 1, a = /3m„ - ^, (2.10) 

where /3 is a free parameter-*^, and rriy = m„/j/2 is the constituent u-quark mass. Neglecting 
the u-baryon production, these N u-quark pairs are assumed to form u-mesons. Typical 
probability distributions of multiplicities are shown in Figs. ^ and ^3 with a Z and Z' 
mediator respectively, for a various values of rrivh- We see the qualitative difference between 
these two regions. For a large values of rriyh, the value of x rrivh is roughly mass of the 
mediator, while for a small values of m^h, N increases approximately logarithmically and 
the value of A^ x m^h is less than the mass of the mediator, with the rest of the energy 
from the mediator decay being carried away as kinetic energy of the light u-mesons. 

In order to determine the momentum of these A'^ v-mesons, we simply adopt the widely 
used Longitudinal Phase Space Approximation (LPSA) which is equivalent to a constant 
distribution of rapidity. The hypothesis of the longitudinal phase space approximation leads 
to a probability distribution which is identical to that of a grand canonical ensemble of non- 
interacting particles. This approximation is used to describe successfully the experimental 
data in e^e~ annihilation, e.g. [28], and is also discussed in many other places, for example 
[29-32]. Since a w-meson i is uniformly distributed in rapidity, its rapidity Yi can be written 
as 

Yi = Z + iiY, (2.11) 

^In this paper, all of the parameters used in the hadronization model are fixed by data from e^e~ 
between -^s =fO and 9f GeV [25]. This has been widely used in the literature. 
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where is a random number between and 1; Z and Y are two arguments which can be 
determined by energy-momentum conservation in the initial vv system: 



N N 



Y,Ei = V^, ^Pu = 0, (2.12) 



i=l i=l 

where Ei and Plz denote the energy and the longitudinal momentum of the ith w-meson. 
They are obtained by 

exp(y,) + exp{-Yi) expjYj) - exp{-Yi} 
Ei = mTi , FLi = mTi , (2-13) 



where the transverse energy rriTi = y nif + P^-, with rrij the mass of the ith f -meson, and 
its transverse momentum Pti obeying the distribution 

N p2 N 

i=l i=l 

a is a free parameter which we again fix with e^e~ annihilation data [25], where a is 
rescaled from the SM value by the ratio Kqcd/^v, where is the confinement scale in 
the valley, assumed to be the T;-hadron mass ra^h in our treatment. Using the simplified 
model described above, we obtain the 4-momenta of the u-mesons. 

2.4 v-hadron decay 

Each ti-hadron produced will subsequently decay to SM leptons or hadrons through the 
diagram shown in Fig. ^o. The branching fractions can be computed from the charges in 
Table ||, with the assumption that 3/4 of the u-hadrons are vectors (which decay democrat- 
ically to heavy and light flavor alike) and 1/4 of the i;-hadrons are pseudoscalars (which, 
like pions, decay predominantly to heavy flavor). The resulting branching to leptons is 
9%, hadrons account for 78% (including ~ 4.7% to r), and missing energy from decays to 
neutrinos about 13%. The lifetimes are computed to be [14]: 

q 1 /10TeV\^ 
r,.^.S~8xl0Vc-'^^^(_) (2.15) 

r,,.4xl0-sec-'-< fl^)'. (2.16) 

(20 GeV)5 \mz'/g' J 

One can see that while the decays are typically prompt, the pseuod-scalar lifetimes 
may be long enough, for lower u-hadron masses, to give rise to a displaced vertex. In 
this case, the displaced vertex is the preferred search method for Hidden Valleys, since 
the backgrounds could be rapidly eliminated. As a result, we focus on the higher mass 
u-hadron case where the lifetimes are sufficiently short that no displaced vertex results, i.e. 
for > 30 GeV. 

In our analyses below, we will decay the ii-mesons to the SM leptons and quarks at 
the parton level. 
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3. Separating Hidden Valleys From Standard Model Backgrounds 



In searching for new physics at hadron colliders, the most important aspect is to identify 
the characteristic and distinctive features of the signal in order to distinguish it from 
standard model processes. Although a concrete prediction of the Hidden Valley signal will 
typically be model-dependent, there are general features that will guide us in the search. 
Generically, the signal events are less jetty, with higher sphericities, cluster masses, and 
typically multiple leptons in the final state, some of which may be isolated. The presence 
of a low mass, narrow resonance in /x pairs from the T;-hadron decay is the most distinctive 
signature in the model we consider. All of these features can be used to reduce the standard 
model background significantly. We turn to quantifying these features step by step. 

As stated in the introduction, for definiteness, we will be studying w-hadrons with 
masses 

m,/,>30 GeV, (3.1) 

presumably determined by the confinement scale, while taking the u-quarks to be essentially 
massless. From Eq. ( p. 15 ) we see that when 't;-hadrons are this heavy, no displaced vertex 



can typically be used to search for u-hadrons. The mass of these hadrons in the absence of 
displaced vertices is typically much larger than standard model hadrons, which have masses 
less than the B meson mass at 5 GeV. To this extent, we have chosen a rather challenging 
f-hadron mass for detection, when there is no displaced vertex for m„/j > 30 GeV. Higher 
mass f-hadrons will be more easily distinguished from the soft QCD background of b, c 
quark hadrons, as we argue in this section. 

Most f-quark and T;-antiquark pairs are produced, as we saw, at the Z or Z' reso- 
nance. In what follows we separate signal on the Z and Z' peaks by the cuts Eqs. ( p.6| ) 
and Eq. ( |2.7D , since the signals are qualitatively different there and they help with the 
background considerations. 

Upon hadronization, an energetic i;-quark gives a broad and massive v-jet on account 
of the higher hadronization scale, with a potentially high multiplicity of u-hadrons in it as 
already shown in Fig. ^. In constructing the events, we take a generous rapidity coverage 
and neglect particles with low transverse momenta. We thus adopt the acceptance for 
particle detection 

\r]\ < 4.9, Pt>3 GeV, (3.2) 

and ignore the charged particles outside this range. Each u-hadron may decay to a pair of 
muons with a branching fraction which is taken to be 4.7%. A pair of muons, even if they 
are fairly soft from the decay of the T;-hadron, can be used as a trigger on the event. We 
take the kinematical acceptance of the two muons as [33] 

Pt(/u) > 10 GeV, \r]{n)\ < 2.5. (3.3) 

In Fig. m, the dashed curve gives the invariant mass distribution of the primary vv after 
the trigger cuts of Eq. (|3.3| ). We see a signal rate reduction of more than an order of 
magnitude near the Z peak, but only a modest reduction near the Z'. This trigger also 
helps considerably with the background separation. We will demonstrate the qualitative 
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Figure 5: pt distribution of the hardest lepton of Hidden Valley events and the QCD backgrounds, 
bh and tt, after the particle acceptance, Eq. (3.2), has been imposed. 



differences between the events from the signal and the SM backgrounds in this section. In 
particular, for the sake of illustration, we will present the QCD background analyses for 
the processes 

pp bbX, ttX, (3.4) 

simulated by PYTHIA [21] with the full QCD initial state radiation, parton showering and 
hadronization. 

The PT distribution of the hardest lepton in the event for the Hidden Valley signal, bb 
and tt events are shown in Fig. ||. In these figures and in what follows, the bb backgrounds 
are sufficiently large that the curves are multiplied by a factor (10~^ in this case) so that 
they can be shown on the same scale as the signal. The presence of the large number of 
soft muons is mainly a result of the collinear behavior of bb along the beam direction and 
the softness of the muons from B decays. 

The pair of leptons from a u-hadron decay form a smaller angle than those from the 
SM Drell-Yan and bb, cc type processes in which the leptons are typically back-to-back in 
the transverse plane. We thus consider a cut for the relative angle in the transverse plane 
between the two hardest leptons in the event (whether they be muons or electrons) 

(pu < vr/2. (3.5) 

This cut will help eliminate most of the lepton pair events from the SM gauge bosons, at 
very little cost to the signal. 

In order to demonstrate the effect of the high confinement scale and hadron mass on 
the shape of the event and on selected observables, such as the thrust, sphericity, and 
cluster mass, we compare in what follows the distributions for a 30 GeV u-hadron mass 
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Figure 6: Differential cross-section as a function of thrust T for rriyh = 30 GcV and the backgrounds 
from bb and it with (a) a Z mediator and (b) a Z' mediator. The sohd hncs indicate the distributions 



before the two muon trigger Eq. ( p.3| ) and the cut on the angle between the leptons, Eq. (3^), are 
imposed, while the dashed lines are after. 

against dominant QCD backgrounds from hh, tt. The collinearity of a jet is often quantified 
in terms of the thrust 

1 = max ' — j — 

_ 22^\Pi\ 

where n is the thrust direction obtained my maximizing the sum. We modify this expression 
in a way that is appropriate for an event in hadronic cohisions, namely to replace all 
momenta of the observed particles with their transverse components 



T = max 



Ei \PTi ■ "I 



(3.6) 



We compare the thrust distributions da/dT in Fig. ^ for a f-hadron signal and bb + X 
background centered on the Z mediator (panel (a)) and tt background on the Z' mediator 
(panel (b)). We show the signal distributions both before (solid curve) and after (dashed 
curve) the two muon cuts Eqs. ( |3.3| ) and ( |3.5| ) are imposed. The backgrounds are shown 
only after these cuts are imposed. We can see that the bb background (and, by extension, 
the background from any low mass hadrons such as D mesons) yields a more jet-like 
structure, with T close to unity. The cut 

T < 0.95 for Z, T < 0.9 for Z' (3.7) 

will very efficiently remove the backgrounds from bb and cc etc. The long tail on the b 
distribution in panel (a) is due to the effects of the initial state radiation. The un-eveness 
in the curves results from a loss of statistics from the stringent muon cut, Eq. ( |3.3D ; this 
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Figure 7: Differential cross-section as a function of sphericity S for ruyh = 30 GeV and the 
backgrounds from bb and tt with a (a) Z mediator and (b) Z' mediator before (sohd) and after 
(dashed) the trigger, Eq. ( |3.3| ) and cut on the angle between the leptons, Eq. (^), are imposed. 



cut reduces an initial event sample of 10 million to just a few tens of events which are 
shown in the dashed bbX curve in the figures. Alternatively, one can define the sphericity 
matrix in the transverse plane 

5(a,/3)^^^, (3.8) 

with a, f3 = x,y. The two-dimensional sphericity (also called circularity) is given in terms 
of the two eigenvalues Qi and Q2 by 

with Qi < Q2 and < S" < 1. We compare the distributions for a u-hadron signal against 
the backgrounds from bb and tt in Fig. ^, again before and after the two muon cuts have been 
imposed. We obtain some marginal improvement for both top and bottom backgrounds, 
where again the un-eveness in the bottom background curve is due to statistics. We make 
a cut 

S>0.1 (3.10) 

for both Z and Z' mediators. 

After obtaining the thrust axis in the transverse plane, we can define two clusters 
around the thrust directions, separated by a longitudinal plane perpendicular to the trans- 
verse thrust axis. The invariant mass of the f-jet cluster will be substantially larger than 
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Figure 8: Differential cross-sections for 30 GeV w-hadron as a function of cluster invariant mass 
with (a) Z mediator and (b) Z' mediator. The hb and tt backgrounds are also shown. All curves 
are after cuts Eqs. (^) and ( |3.5[ ) have been imposed. 

that of a standard model jet. We sum all observable particles in a cluster 

The cluster mass will typically be half the mediator mass for a fairly spherical event, becom- 
ing smaller as the event becomes increasingly jetty. We show the cluster mass distributions 
in Fig. ^ for Z mediator (panel a) and Z' mediator (panel b). The tt and 66 backgrounds 



after cuts of Eqs. (3^) and (|3.5| ) are also shown. In contrast to a Hidden Valley event 
where the cluster mass is of the order of a half of the mediator mass, the cluster invariant 
mass for the heavy quarks should be peaked near Mduster ~ f^t- Thus we can further 
remove the heavy quark backgrounds with a simple cut 

M,.... > I (3.12) 
I 200 GeV at Z'. ^ ' 

We would like to point out another important feature for the signal events. The two clusters 
should have roughly equal cluster mass up to the fluctuations of the v and v hadronization 
and the decay of the v mesons. It can be seen from Fig. ^ that the mass spread is only about 
30ibl0 GeV at the Z pole and 250ib50 GeV at the Z' . For higher order backgrounds on the 
other hand, such as /i^^^+jets, cc, 66-|-jets, the event shapes would be rather asymmetric 
on both sides. 

In order to develop an intuitive picture of a typical event satisfying the basic cuts, we 
show the lego plots for two typical events satisfying the basic cuts Eqs. ( |3.3D and ( |3.5D in 
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Figure 9: Lego plots in the rj — (j> plane for typical events satisfying the basic cuts Eqs. (3^) and 



(3.5) for (a) a Z mediator and (b) a Z' mediator. The color codes indicate the hadronic energy 
(red), electromagnetic (blue), and muons (green). Hadronic energy is at the parton level only (no 
showering). 



Fig. ^ for one mediator each. The height indicates the energy scale for hadronic energy 
(red), electromagnetic (blue), and muons (green). As expected, these events are much 
fatter than QCD events, with high multiplicities, and often many leptons in the final state. 
The plots for the signal are at the lepton and quark level, without QCD hadronization. 
In a more realistic situation, the energy spread will be even broader. To contrast Hidden 
Valley events with the backgrounds in the typical case, bb, tt we consider here, we also 
show two typical events satisfying the same basic cuts in Fig. We see clearly the more 
isolated jetty structure. The background events are generated by PYTHIA [21] with the 
full QCD initial state radiation, parton showering and hadronization. 

With this set of judicious cuts, we can now proceed to the decisive step, to reconstruct 
the u-hadron mass. This reconstruction can provide ultimate power to distinguish the 
signal from the backgrounds. We propose to reconstruct the 't;-hadron resonance from the 
leptonic decay products, in particular a pair of muons. Quite often, there are more than 2 
muons in an event (n^ ^u^'s and n_ /^~'s). We form all the n+n_ combinations of /x^/i^. 
Among all possible opposite sign pairs in the event, if there are k pairs with the same 
invariant mass (±1 GeV), we weight each pair with l/fc, throwing the pairs away that do 
not match. In this way we are able to select the correct muons to pair together. In the case 
that we are missing a muon (from the detector acceptance) and none of the pairs match, we 
compute all possible pairs and weight it by l/(n_|_?7-_). This is a combinatorial background 



that reduces the signal. This is shown in Fig. 11, for both the Z (panel (a)) and Z' (panel 
(b)) mediators. With the excellent momentum resolution at these low muon energies, we 
take 500 MeV bins, consistent with experimental resolution at those energies [33], and 
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Figure 10: Lego plots in the rj — cf) plane for typical (a) bb events (on the Z-mediator) and (b) tt 



events (on the Z'-niediator) satisfying the basic cuts Eqs. (O) and (1_5). The color codes indicate 



the hadronic energy (red), electromagnetic (blue), and muons (green). 

show the number of events in each bin predicted for 100 fb~^ (on Z peak) and 10 fb~^ (on 
Z' peak) of data. We also show the results from bb, tt backgrounds in these two panels 
respectively, which demonstrate the final effectiveness to separate the backgrounds. The 
solid lines (for both signal and background) are with the 2fj, trigger and cpn cuts only; the 
dashed lines (both signal and background) include all cuts in Table 2 through M duster- It 



is noted from Fig. 11(a) that the background muon pairs (solid and dashed) are all from 
the 6-quark cascade decays. This implies that the higher order contributions of high pT 
bb events may not have been adequately generated in the PYTHIA simulation. We thus 
perform a parton level calculation for bb plus a high pT jet with bb — > ji^ iJi~X. Its mn 



distributions are shown in Fig. 11 for both Z and Z' mediators, with only the basic trigger 
implemented. The shape variable cuts {S, T and Mduster) cannot be properly implemented 
on this background since it is computed at the parton level, although such cuts may be 
very effective in reducing it. We can impose a final resonant mass cut to estimate the signal 
significance 

rriyh — 1 GeV < mu < m^h + 1 GeV, (3.13) 

where mu is the invariant mass of the lepton pair selected for the event. Often there is 
more than one muon resonance in the event which can even further strengthen the signal 
identification. 

We have now finished laying out the cuts which will distinguish the Hidden Valley 
Model from the Standard Model events. We summarize these cuts, laid out step- by-step 
in Eqs. (U), (jsj), (|3), (|3lo| ), (|1|), (|1|), in Table |, where we show their effect on 



the signal rate. The entries in the last row show the tt background with the consecutive 
cuts corresponding to the signal of a Z' mediator. This background becomes much smaller 
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Figure 11: The mass reconstruction of tlie u-liadron via the lepton pairs in the decay, (a) for the 
Z mediator and (b) for the Z' mediator. The bin size is 500 MeV, consistent with experimental 
resolution at those energies. The solid lines (for both signal and background) are with the 2[i trigger 
and 4>ii cuts only; the dashed lines (both signal and background) include all cuts in Table 2 through 
Mciuster- The background curves for 66+jets are before the cuts on S,T are imposed. 



than the expected signal. We have not included the backgrounds for the case of a Z 
mediator, since there are subleading, but possibly significant processes like bb+ jets that 
have not been adequately simulated in PYTHIA with showering. Further work is required 
to show quantitatively and conclusively that the signal can be dug out of the background 
on the Z peak. By contrast, we have shown that the signal should be easily separable from 
background on the Z' peak. 

The cuts that we have specified to this point show the main features of the Hidden 
Valley signal which will separate it from background. There is an additional feature, 
however, which may also be useful for the signal separation should it be necessary. When 
the ti-hadrons decay to pairs of standard model particles via the mediators, the separation 
between the decay products (a pair of leptons or a pair of quarks) will be larger than 
is typical of standard model hadron decays, since the angular size of a typical decay is 
m^uhf^Efi, where Eh is the energy of the t;-hadron. As a result of the more spherical jet 
structure of Hidden Valley events and the larger separation between decay products, an 
isolated lepton can be a feature to select for Hidden Valley events. This is demonstrated in 
Fig. |l2|, where Ai2^°^ is the separation between the most isolated lepton (whether muon 
or electron) and its nearest (non-leptonic) neighbor. An additional cut on the leptons may 
then be designed so that 

PTW>6GeV, |r/(£)|<2.5, Afij^°i > 0.3. (3.14) 

This will substantially remove the SM backgrounds, especially those from heavy quarks, 
while necessarily reduce the signal rate as well, depending on the number of isolated leptons 
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Figure 12: Differential cross-section as a function of the separation for tlie most isolated lepton 
ARl^°l^ for iTLyh = 30 GeV with (a) a Z mediator and (b) a Z' mediator. 
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Table 2: Signal cross-sections (in fb) near the peaks of the Z and Z' mediators with various 
kinematical cuts (defined in the marked equations in the text). Cuts up to rriyh are implemented 
cumulatively, while the cuts to tlie right of the double line are implemented independently of each 
other, but in tandem with the basic cuts (those cuts through Mduster)- We also show the effect of 
the cuts on the dominant background tiX for Z' mediator. 



required, as can be seen from the last three columns in Table Q before imposing the niyh 
cut. 

We hope that by now we have demonstrated enough characteristics for the Hidden 
Valley signal and have developed the judicious cuts to separate it from the SM backgrounds. 
We have treated in detail the backgrounds which are the most dangerous. We summarize 
now the potential backgrounds, with reference to how the cuts we have defined in this 
section have helped to remove them. Of course, all of these backgrounds are continuum in 
m^^ distribution, above which a new narrow resonance may easily stand out. 

• Drell-Yan background, that is, Z and 7* plus jets, is of order 1 nb, and 

is reduced on the Z' by the requirement of (^]^) to the order of 7 fb. To obtain 
two hemispheres of large cluster mass requires of order three additional hard jets. 
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reducing the background below problematic levels. W plus heavy flavor can similarly 
produce dileptons, but is similarly reduced. 

• Gauge-boson pair background: Another source of di-muon events may come from 
the gauge boson pair production W~^W~ ,WZ, ZZ — > ^~ +X. The production 
rates for these processes start at the order of 100 pb, and are reduced to the order 
of 500 fb on the Z' . Dimuon branching fractions reduce this by almost two orders of 
magnitude, and to obtain high cluster mass on both sides requires additional jets. 

• cc, bh backgrounds: A more difficult background results from the QCD strong pro- 
duction of heavy quarks cc and hb that can lead to a large cross section of the order 
1 /xb at pt ~ 30 GeV. With their semi-leptonic decays to fi~^n~ with approximately 
10% branching fraction and the basic triggering cuts, the rate will again be at a level 
of 1 nb. Similar considerations as those to eliminate the DY background, along with 
an additional requirement that the system mass is near either the Z or Z' mediator, 
two massive clusters, and anti-tagging for c, b will reduce this background by many 
orders of magnitude. We have shown by PYTHIA simulations that the leading order 
in Os hb background process could be effectively separated from the signal, as seen in 

Fig. (InD. 

Additional radiation of jets or multiple heavy quarks may make the event topology 
more spherical. However, the qualitative features, such as the jetty structure, non- 
equal cluster masses, lepton isolation etc., are expected to provide additional handles 
to discriminate those higher order backgrounds. For illustration, we calculated the 
process of hb plus a high px jet with bh fi^fi^X in Fig. (0). This background 
is under control at the Z' peak, but sizable at the Z peak. Nevertheless, we expect 
that it may be sufficiently suppressed by the sphericity, thrust, cluster mass, and in 
particular the lepton isolation, though simulations beyond the scope of this paper 
would be required to confirm this. 

• tt background: The top-quark pair production at the LHC has a total cross section 
of 800 pb, and is reduced on the Z' to about 20 pb. Di-muon decays of the top 
quarks cost another factor of 80, but typically the angle between the two muons is 
very large, since the top quarks are boosted. The most likely way to obtain two 
muons passing our cuts is from t — > Wb in which both W and h produce a muon; 
since this can happen for either i or f it costs a factor of about 40. With the cluster 
mass requirement, our study has shown that the dimuon invariant-mass continuum 
background is small. Additional jets can push the event above the cluster mass 
requirements, but the signal is so far above background that it will still easily stand 
out. 

Since we have already seen that our signal stands out far above the tt background, we 
are confident that inclusion of a larger set of continuum backgrounds and of detector 
backgrounds will not change our basic conclusion. If it does, then a lepton isolation cut 
will eliminate the largest backgrounds, without costing too much of the signal, as is clear 
from Table 2. 
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We summarize the total cross-section for the signal in Table ^ including all the cuts 
and leptonic selections step by step. Based on the above arguments on the background 
rates and their possible suppressions, we are led to believe that the SM backgrounds are 
effectively removed. The signal rate at the Z' resonance is higher than that at the Z by 



two orders of magnitude with out parameter choices of Eq. (2.3). We can thus estimate the 
sensitivity for the signal observation for different values of the parameters. Assuming an 
integrated luminosity of 10 fb~^ and take 4 background events as in Table H this requires 
15 signal events to reach a 99% CL observation with Poisson statistics. We can thus reach 
a coupling strength 

g' ~ 1.7 X 10~2 for mz' = 1 TeV. (3.15) 
Conversely, if we fix the coupling, we obtain a 99% CL sensitivity on the mass 

mz' ~ 4.6 TeV, for g' ~ 1/7, (3.16) 

after including the fall of the valence quark luminosity at high momentum fraction x values, 
assumed to go like (1 — xi)^(l — ^2)*^. 

We reiterate that we have not performed detailed simulations including the realistic 
experimental environment and detector effects, which would be necessary to draw a more 
definitive conclusion. 



4. Summary and Conclusions 

In this paper, we have studied the feasibility of detecting hidden sectors with low mass 
bound-states at hadron colliders. We assumed a generic gauge boson Z' at the TeV scale 
as the mediator (that also mixes with the SM Z) to induce the interactions with the SM 
sector. Most of the Hidden Valley signal results from a Z or Z' produced on-peak which 
then decays to Hidden Valley f-quarks. The t;-quarks subsequently confine into u-hadrons, 
which then decay to SM pairs of quarks or leptons. We chose the characteristic case with 
a pair of muons in its decay. We took the confinement scale to be m^fi > 30 GeV, below 
which the f-hadrons may be long-lived and thus would lead to different signatures. We 
demonstrated the characteristic features of the signal with a Monte Carlo simulation of the 
hadronization, and contrasted them with the SM background expectations. We showed 
that accounting for the constraint from the Z pole physics at LEP, searches at Tevatron 
may be difficult in this particular class of Hidden Valleys. We thus concentrated on the 
searches at the LHC. We found that for the signal events 

• there may be multiple hard muons {pT > 10 GeV), more widely separated, which 
can be used as triggers, as seen in Figs. ^ and 12; 



the events are more spherical and less thrusty, as shown in Figs. ^, |^ and 

the invariant mass of the f-jet is high: u-jets are very "fat", and a signal event has a 
pair of back-to-back clusters of nearly equal mass, as shown in Fig. & 
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muon pairs from the spin-one i;-meson decay can reconstruct niyh as a narrow reso- 
nance, and provide the most effective kinematical variable to identify the new signal, 



as shown in Fig. 11. Often there is more than one such resonance per event. 



The signal rate near Z' is substantially larger than that at Z, especially after imposing 
acceptance cuts on the leptons. We have designed judicious cuts to remove or reduce the 
SM background, the most difficult of which are bb at high pT, but relatively low center of 
mass energy (near the Z peak). The signal on the Z peak is much more challenging than 
on the Z' peak, due to the larger b quark backgrounds. By contrast, we have shown that 
the tt backgrounds on the Z' mediator can be easily and efficiently removed. On the Z', 
we found large coverage for the Hidden Valley parameter space {g' , mz') at the LHC as 
shown in Eqs. ( p. 15 ) and ( 3.16| ). 



We emphasize that although we studied a particular class of Hidden Valley models, 
the implications of our efforts are much more general. Light spin-one resonances that can 
decay to dilepton pairs, and that are produced only in rare decays of the Z or in decays 
of a new particle such as a Higgs or Z', may easily arise in a very wide array of Hidden 
Valley models. They may also arise in yet other classes of models. A search for such 
resonances cannot succeed in a fully inclusive dimuon channel. Instead, one must select 
events on the basis of variables such as total visible energy and event-shape observables; 
other examples might include missing energy, numbers of b tags, etc. These variables cut 
away most standard model processes, and with the events that remain, narrow dilepton 
resonances stand out easily, even for signals with small cross-sections. We view this as a 
general lesson for LHC, and even Tevatron, search strategies. 

We have studied the basic features of models with light quarks, confinement scale above 
~ 10 GeV, a light spin-one resonance, and a gauge mediator. We note that there are other 
Hidden Valleys which remain to be explored. The phenomenology depends on the matter 
content of the hidden sector, the size of the mass gap (confinement scale), and the nature 
of the mediator, each of which may have unique implications for general collider searches 
for Hidden Valleys, as well as for astrophysics and cosmology. With these tools to study 
low mass hidden sectors at hadron colliders, we may uncover a richer sector beyond the 
mostly discussed new physics at the TeV scale. 
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